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We consider a 4D model for photon production induced by a refractive index perturbation in a 
dielectric medium. We show that, in this model, we can infer the presence of a Hawking type effect. 
This prediction shows up both in the analogue Hawking framework, which is implemented in the 
pulse frame and relies on the peculiar properties of the effective geometry in which quantum fields 
propagate, as well as in the laboratory frame, through standard quantum field theory calculations. 
Effects of optical dispersion are also taken into account, and are shown to provide a limited energy 
bandwidth for the emission of Hawking radiation. 



I. INTRODUCTION 

One of the most intriguing predictions of quantum 
fields in curved spacetime geometries is the production 
of Hawking radiation. In 1974 S. Hawking predicted 
that black holes emit particles with a thermal spectrum. 
Therefore a black hole will evaporate, shedding energy 
under the form of a blackbody emission [I}j3]. However, 
it turns out that for a typical stellar mass black hole the 
temperature of this radiation is so low (~10 nK) that it 
has no hope of being directly detected. Notwithstand- 
ing, Hawking radiation does not actually require a true 
gravitational event horizon but, rather, it is essentially a 
kinematical effect, i.e. it requires some basic "kinemati- 
cal" ingredients but no specific underlying dynamics (see 
e.g. [IHS]). What suffices is the formation of a trapping 
horizon in a curved spacetime metric of any kind, and the 
analysis of the behavior of any quantum field therein (see 
also the discussion in [6]). The quanta of this field will 
then be excited according to the prediction of Hawking. 
On this basis, a number of analogue systems have been 
proposed, for the first time by W. Unruh [7] and later 
by other researchers (see e.g. [8] and references therein), 
that aim at reproducing some aspects of the kinematics 
of gravitational systems. Most of these analogies rely on 
acoustic perturbations and on the realization of so-called 
dumb holes: a liquid or gas medium is made to flow faster 
than the velocity of acoustic waves in the same medium 
so that at the transition point between sub and super- 
sonic flow, a trapping horizon is formed that may be tra- 
versed by the acoustic quanta, viz. phonons, only in one 
direction [5]. Unfortunately the phonon blackbody spec- 
trum is expected to still have very low temperatures, thus 
eluding direct detection (see e.g. [S]). A parallel line of 
investigation involves effective geometry for light, which 
has been introduced by Gordon [TU] and extended also to 
nonlinear electrodynamics; black hole metrics have been 
introduced and the possibility to perform experiments 
involving Hawking analogue radiation has been explored 
[TlTf2"5] . Some time ago Philbin et al. proposed an optical 



analogue in which a soliton with intensity /, propagating 
in an optical fibre, generates through the nonlinear Kerr 
effect a refractive index perturbation (RIP), Sn = ri2l, 
where 112 is the Kerr index 24J. The same mechanism 
has also been generalized to a full 4D geometry by Fac- 
cio et al. [25] and has led to the first experimental ob- 
servation of quanta emitted from an analogue horizon 
|26| . The RIP modifies the spacetime geometry as seen 
by co-propagating light rays and, similarly to the acous- 
tic analogy, if the RIP is locally superluminal, i.e. if it 
locally travels faster than the phase velocity of light in 
the medium, an horizon is formed and Hawking radiation 
is to be expected. 

Here we take into consideration the Hawking effect in 
dielectric black holes induced by a RIP which propagates 
with constant velocity v. We first show that, under- 
suitable conditions, we can re-map the original Maxwell 
equations for nonlinear optics into a geometrical descrip- 
tion, in analogy to what is done in the case of acoustic 
perturbations in condensed matter. Then we provide a 
model in which the presence of the Hawking radiation 
can be inferred even without recourse to the analogous 
model characterized by a curved spacetime geometry. 
Indeed, Hawking radiation is a new phenomenon for 
nonlinear optics, never foreseen before, which could 
legitimately meet a sceptical attitude by the nonlinear 
optics community, because of a missing description of the 
phenomenon by means of the standard tools of quantum 
electrodynamics. As a consequence, it is important to 
recast the analogue picture by using "standard" tools 
of quantum field theory (without any reference to the 
geometrical picture). 

The peculiar feature which distinguishes the dielectric 
Hawking effect from the traditional black hole one is the 
presence of optical dispersion, which gives rise to relevant 
physical consequences on the quantum phenomenon of 
particle creation. Therefore, we take dispersion into 
account. We first introduce a 2D reduction of the model 
in the presence of dispersion, and see how the dispersion 
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relation is affected by the frequency-dependence of 
the refractive index. Then we discuss both phase 
velocity horizons and group velocity horizons, which can 
both play a relevant role in the physical situation at 
hand; their qualitative difference and the possibility to 
discriminate between them in experiments (numerical 
and/or laboratory-based) is considered. 



Section |V] is devoted to the conclusions. 

Finally, we have added three appendices. The first es- 
tablishes a correspondence between our black hole metric 
and the acoustic one. The second one relates the temper- 
ature to the standard conical singularity of the Euclidean 
version of the metric. The third establishes the relevant 
analytic properties of the Bogoliubov coefficients. 



The structure of the paper is as follows. In section |TT] 
we introduce our idealized model of a RIP propagating, 
in a nonlinear Kerr medium, with a constant velocity v 
in the x direction and infinitely extended in the trans- 
verse y and z directions. In the eikonal approximation, 
the model is embodied in a suitable wave equation for 
the generic component of the electric field propagating 
in the nonlinear medium affected by the RIP. We de- 
scribe this propagation as taking place in an analogue 
spacetime metric, written in the pulse frame, where the 
metric is static and displays two horizons x+ and X- 
for propagating photons, analogues to a black hole and 
a white hole horizons respectively. By assigning to the 
black hole horizon a surface gravity in a standard way, 
we calculate the evaporation temperature T + of the RIP 
in the latter frame, which turns out to be proportional to 
the absolute value of the derivative of the RIP's profile, 
evaluated at x + . We also show that T + is conformally 
invariant, as expected for the consistency of the model. 

In section [nil we tackle the model from a different per- 
spective, namely within the framework of quantum field 
theory. We construct a complete set of positive frequency 
solutions of the wave equation in the lab frame, in 
terms of which we quantize the field in the standard way. 
Then, we compare the ^k, with the corresponding plane 
wave solutions of the wave equation in the absence of the 
RIP. This enables us to calculate the Bogoliubov coeffi- 
cients, in terms of which we express the expectation value 
of the number operator of the outgoing quanta in the in 
vacuum state. As expected, this average is a thermal-like 
distribution displaying the typical cos 9 temperature pro- 
file in the lab frame, where 9 is the emission angle with 
respect to the direction of motion of the RIP. We then 
display the transformation law which relates the temper- 
atures in the pulse and in the lab frame, respectively, by 
means of the usual Dopplcr formula. 

The results worked out in the preceding sections apply 
in the approximation in which dispersion is neglected. In 
section IV we study the modifications we expect when 



dispersion is taken into account. In particular, we show 
that, whereas in the dispersionless case one should expect 
to observe all the blackbody spectrum, in the presence 
of dispersion only a limited spectral region, which de- 
pends on v, will be excited, and the blackbody spectrum 
shape will not be discernible. Moreover, one can intro- 
duce two different concepts of horizon: phase velocity 
horizons and group velocity horizons, which are shown 
to be involved with different spectral bandwidths, and 
with very different qualitative behavior regarding their 
action on photons. 



II. STATIC DIELECTRIC BLACK HOLE 

We consider a model-equation which is derived from 
nonlinear electrodynamics (with \2 = and X3 0) in 
the eikonal approximation for a perturbation of a full- 
nonlinear electric field propagating in a nonlinear Kerr 
medium. In order to make the forthcoming analysis as 
simple as possible we replace the electric field with a 
scalar field $ (cf. e.g. Schwinger's analysis of sonolu- 
minescence |28j). for which we obtain the wave equation 



n 2 (xi - vti) 2 



oft - df$ = o. 



(i) 



Coordinates in the lab frame are denoted by ti,xi,y,z 
(the labels of y, z are omitted because they will not be 
involved in the boost relating the lab frame with the pulse 
frame). Here, n{x\ — vti) is the refractive index, which 
accounts for the propagating RIP in the dielectric. The 
latter can be obtained by means of an intense laser pulse 
in a nonlinear dielectric medium (Kerr effect). In our 
model the RIP does not depend on the transverse coor- 
dinates, namely it is infinitely extended in the transverse 
dimensions. This approximation is necessary in order 
to carry out the calculations below, which allow us to 
draw a tight analogy between the Hawking geometrical 
description and the quantum field theoretical treatment. 
Such an infinitely extended RIP is clearly an idealization, 
which cannot be obtained in an actual experiment. How- 
ever, relatively flat, e.g. super-Gaussian-like pulses may 
be produced that would, at least locally, fall within the 
approximations adopted here. 

Eq. arises also in the eikonal approximation for a 
scalar field in the metric 



ds 2 



n 2 (xi — vti) 



dxf 



dy — dz 



(2) 



To carry out the analysis in the context of the analogue 
metric, we pass from the laboratory frame to the pulse 
frame (refractive index perturbation frame) by means of 
a boost: t = j(ti — jsXi), x — j(xi — vti), and we obtain 
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FIG. 1: (in color online) Scheme of the RIP geometry. x+ and 
x- indicate the black hole and white hole horizon positions, 
respectively. 



We assume 



n(x) = no + Sn = no + t]I(x), 



(4) 



where r\ is meant to be positive, with 77 <§C 1 due to the 
actual smallness of the Kerr index (a negative r\ could be 
easily taken into account); I denotes the normalized in- 
tensity of the pulse, is taken to be a C°° function, rapidly 
decaying at infinity and with a single maximum at x = 0, 
of height 1. A scheme of the RIP is shown in Fig. [I] The 
form of I implies that both dt =: £ and d$ are Killing 
vectors for the given metric. Then the surface goo = is 
lightlike and corresponds to an event horizon. It is also 
possible to verify that the Frobenius integrability condi- 
tions are trivially satisfied, so that there exists a coor- 
dinate transformation carrying the metric into a static 
form (see e.g. [30 ). Even if these coordinates arc sin- 
gular, we carry out the relative transformation because, 
on one hand, it allows a more straightforward comparison 
with the well-known Schwarzschild case, and on the other 
hand, it allows a direct computation of the greybody coef- 
ficient (cf. sect. HID |. To implement this transformation 



consider the following coordinate change: 



where 



dt = dr — a(x)dx, 

, s 5oi 0*0 
a{x) = — 

goo{x) 

Then the metric takes the static form 

1 



(5) 



(6) 



ds 2 = 



n 2 (x) 



g TT (x)dr ~ 



9tt{x) 



dx 2 -dy 2 -dz 2 , (7) 



where 



g TT (x) := 7 2 + n(x)^j (l - n(x)^j 



(8) 



There is a remarkable resemblance of the r, x-part of the 
metric with a standard static spherically symmetric met- 
ric in general relativity in the so-called Schwarzschild 
gauge, aside from the important difference represented 
by the factor ^ replacing c . The horizons are deter- 
mined by the condition g TT = 0, i.e. by 



which, by Eq. Q, becomes no + r]I(x) = c/v. Due to the 
specific shape assumed for / there will be two horizons: 
one, denoted by x + , located on the rising edge of the 
RIP (i.e. dn/dx\ x+ < 0) and one, denoted by X—, on the 
falling edge (i.e. dn/dx\ x _ > 0). Since < I < 1, the 
condition for the occurrence of the event horizons is 



1 



n + r) 



v 1 

< - < — 
c 



"0 



(10) 



The external region corresponds to 1 < i_ and to 
x > x+. The leading horizon x — x+ is a black hole 
horizon, whereas the trailing one x = x_ is a white hole 
horizon. This can be understood as follows. Consider the 
front observer, i.e. the observer toward which the pulse 
is moving. Any photon starting from the region behind 
the front zero x+ cannot reach this observer before the 
latter is reached by the pulse. Indeed, the pulse is mov- 
ing superluminally relative to any photon travelling in 
the region between the two zeros x_ and x+. Therefore, 
the "leading" zero x+ corresponds to a black hole hori- 
zon. Analogously, the "trailing" zero X- corresponds to 
a white hole horizon. See also Appendix [X] for further 
details. 



A. Surface gravity and temperature 

We can formally assign a temperature to the "black 
hole" horizon z + by defining it in the usual way in terms 
of a "surface gravity" k + associated to the latter. We 
have 

1 
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dx J 



(11) 



X=X + 



or [ST] 



2 2 

«+ = 7 v 



dn 
dx 



(*+) 



(12) 



Then, the temperature is given by the familiar formula 
which has been also derived in 1241: 



2 2 
7 v 



27rfc;,c 



dn 
dx 



0*+) 



(13) 



We can arrive at formula ( 12 ) by several different meth- 
ods. For the sake of completeness, in Appendix [B 



we 



employ a different method to derive the expression ( 13 ) 
for the temperature, and its invariance with respect to 
the conformal factor. The latter aspect is well-known in 
General Relativity [3T] , and also in the frame of acoustic 
black holes [32]. 
Since 



1 — n(x) - 



0. 



(9) 



n(x + ) — n + kr], 



(14) 
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with k = I(x+) G (0,1), we can rewrite equations (12) 
and (fl3| respectively as follows: 



(no + ki]) 2 - 
he 1 
27rfcb (no + krjy 



dn 
dx 



(*+) 



dn 
dx 



(*+) 



(15) 
(16) 



B. Gaussian pulse 

As an example of a refractive index perturbation r)I(x) 
consider a Gaussian normalized intensity of the pulse, 



I(x) = exp — 



x 2 \ 
2a 2 ) 



(17) 



and r\ a small parameter. These choice reflects a typical 
situation in experimental optics. Then the horizons given 
by Eq. Q ) occur for 



x± = ±<7<t 







/-21og 


(- -n ) - 







(18) 



and X- = —x+. The surface gravity and the temperature 
are respectively given by 



(19) 



and by 



|24) ) . This is a typical e.g. of spontaneous laser pulse 
dynamics in transparent media with third order nonlin- 
earity. An analogous behavior is expected for the refrac- 
tive index. As a consequence, in our model it is necessary 
to introduce, beyond the scale 2a, which, roughly speak- 
ing, represents the overall spatial extension of the RIP, at 
least a further scale that corresponds to the 'thickness' of 
the region where the refractive index undergoes its most 
rapid variation. Then we can adopt a profile which is 
analytically described as follows: let us define 



H(x) := 1 + tanh 



■>wh 



tanh 



(22) 



where S w h, S^h > are length scales describing the 'thik- 
ness' of the region over which a rapid variation of the 
refractive index occurs, and 2a > corresponds to the 
overall spatial extension of the RIP. Compare also the 
choice of the velocity profile for acoustic geometries cho- 
sen in |33j . Then for the refractive index profile we choose 



n(x) =no + r] 



H(x) 



max x H(x) ' 



(23) 



It is evident that asymptotically n converges to hq and 
the RIP is correctly normalized. This model can be also 
useful in view of exploring the possibility to obtain a 
bh-wh laser [33] . 

From a physical point of view, we point out that, in the 
presence of a shock front typical values of S w h,T] may 
be order of 6 w h ~ 10 -6 m and rj ~ 5 • 10 -3 respectively, 
leading to T + ~ 1 K. 



T, = 



2irkhC 



he 1 



(20) 



Typical values of the parameters are: a ~ 10 5 m, r\ 

10- 



3 , n ~ 1.45. Then, for k — §, say, we get 



2 ■ 10~ 2 K, 



(21) 



which is greater than in standard cases of sonic black 
holes. This temperature value is not to be meant as char- 
acteristic, somewhat larger values may be obtained. See 
the following subsection. 

Note that this is the temperature in the frame of the 
pulse, with which the geometry is associated. A further 
transformation is needed in order to recover the temper- 
ature in the laboratory frame (see the following section). 
Note also that the temperature is proportional to the 
derivative with respect to x of the refractive index, eval- 
uated at the black hole horizon. This is in agreement 
with the results of Ref. [Ml. 



C. Shockwave model 

Let us assume that the rear part (trailing edge) of the 
filament is characterized by a Shockwave profile (see e.g. 



III. QUANTUM FIELD THEORETICAL 
TREATMENT: HAWKING RADIATION IN THE 
LABORATORY FRAME 

The deduction of Hawking radiation emission by a 
(non-extremal) black hole can also be carried out in a dy- 
namical situation where the effects of the onset of a black 
hole horizon on quantum field theory modes is taken into 
account. See e.g. the seminal paper by Hawking [5] and 
also HTj: where the equivalence between the static and 
the dynamical picture for Hawking radiation is consid- 
ered. In the following section, we develop a dynamical 
picture for particle production, with the aim of not re- 
lying only on the analogue gravity picture but also ob- 
taining the phenomenon in a non geometrical setting by 
means of standard tools of quantum field theory. 
Consider a massless scalar field propagating in a dielec- 
tric medium and satisfying Eq. ([l]). It is not difficult to 
show that this equation arises in the eikonal approxima- 
tion for the components of the electric field perturbation 
in a nonlinear Kerr medium |12l 115] . n(xi — vti) is the 
refractive index, which accounts for the propagating re- 
fractive index perturbation in the dielectric. 

In order to develop the model we adopt the following 
strategy: a) we look first for a complete set of solutions 
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of equation ([lj; b) then we employ these solutions to 
perform a comparison between an initial situation, con- 
sisting of an unperturbed dielectric without a laser pulse, 
with a uniform and constant refractive index hq, and a 
final situation where a laser pulse inducing a superlumi- 
nal RIP is instead present. This scenario is analogous to 
the original situation envisaged and treated by Hawking 
[5] in which one considers an initial spherically symmetric 
astrophysical object, with no particles present at infinity, 
followed by a collapsing phase leading eventually to an 
evaporating Schwarzschild black hole. In analogy with 
the original treatment by Hawking, our aim is to evalu- 
ate the mean value of the number of quanta (calculated 
by means of final creation and annihilation operators) on 
the initial vacuum state. 



A. Out states 

Recast equation ([!]) in terms of the new retarded and 
advanced variables, respectively u = xi — vti and w — 
xi + vtf. 



n 2 {u)v 2 



(d 2 u <s> + dl<s>-2d u d w <s>) 



- + €® + - - flg$ = 0, 

and look for monochromatic solutions of the form 



(24) 



(25) 



Then, the ansatz of Eq. (25) leads to 



A''(u) + 2ik w C l + n l { ' u) f 2 A'(u)- 

_ ( ,2 , k l + k z 



A{u) = 0. (26) 



Assuming n(u) to be analytic (see e.g. (17)) we see that 



the coefficients of A'(u) and of A(u) have a first order 
pole at the roots of 



1 - n(u)- = 0. 



(27) 



Then, assuming n to be of the form Q, the condition 
for the occurrence of the event horizons is given by (10). 



Since x — 71*, the black hole and white hole horizons, 
when they exist, are located respectively at 



u± 



1 

7' 



-x±, 



(28) 



and we have a second order linear differential equation 
with Fuchsian singular points at u = u± (with u_ < 
The general solution of Eq. ( 24 ) can be determined, in the 



neighbourhood of the singular points u±, by the stan- 
dard methods of integration by series. We perform this 
calculation in Appendix [Cl Here, instead, since we are 



interested in the asymptotic behaviour of the solutions 
as u — >• u± and as u — > ±00, the best strategy is to ap- 
peal to the WKB approximation. Namely, we write $ as 
in Eq.p5l with 



A(u) — exp I i / k u {u')du' 



(29) 

Then, we find the following dispersion relation: 

7l 2 V 2 

s (k u ~~ k w ) — (k u + k w ) — kj^ = 0, (30) 
c l 

where k\ = k 2 + k 2 . By solving for k u as a function 
of k w , kj_ one arrives at a second degree equation whose 
solutions arc of the form: 



1 



-± 



±2- 



k\ 1 



kl 4i 



The constraint 



1,2 

n ,,, 



< 



1 



(31) 



(32) 



must hold in order that propagating solutions are avail- 
able. This seemingly would imply the existence of a lim- 
iting angle for the emission, but one has to take into 
consideration that, in approaching the horizon, the right 
hand side of the latest equationn becomes infinite, leav- 
ing room for no real limitation on the emission angle. 

The root fc+ is singular at the horizons, whereas k~ is 
regular there. We focus on the behaviour of fc+ in the 
neighbourhood of u + . Near u + the dependence on kj^ is 
washed out, in agreement with the analysis in [5 , and we 
have 



2&71 



Then, 



A(u) ~ exp 



vn'(u + ) 



log(u - U+) 



(33) 



(34) 



The logarithmic divergence in the phase as u \. u + is en- 
tirely analogous to the one experienced by the covariant 
wave equation at the horizon of a Schwarzschild black 
hole [5J [3] and represents the typical behaviour of out- 
going modes approaching the horizon 5J. It is worth to 
remark that this divergence appears here without any ref- 
erence to geometry. In addition, as shown in Appendix [C| 
it emerges naturally as an exact asymptotic (as u \. u + ) 
consequence of Eq. ( 24 ) . On the other hand, the same 



asymptotic behaviour can be inferred also in the context 
of the geometric approach. Indeed, such behaviour is 
characteristic also of the monochromatic solutions of the 
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covariant wave equation in the metric though the 
latter equation is not exactly solvable. 

Regarding the dependence on k± , we note that it does 
not disappear far away from the horizon, but this turns 
out not to be relevant for the Hawking phenomenon [5]. 
Also note that k u has the opposite sign of k w , and \k u \ > 
I k w | . 

We denote by F + the solutions associated to the root 
k+ . Their asymptotic form (at large distances from the 
pulse, i.e. for u 3> u + ) is 



^ £,i(k u J rk lu )xi-\-ik y y-)-ik z z — ivti(k u — k 1L <) 



This can be written as 



F + = F+ ~ exp (iki ■ x ; - iuiU) . 



where k ; = {k x i,k y ,k z ), with 



uji = v(k u - k w ), 



so that 



k u 
u 



(35) 



(36) 



(37) 



(38) 
(39) 



(40) 
(41) 



We look for solutions which are asymptotically of posi- 
tive frequency and outgoing with respect to the dielectric 
perturbation, which is right moving. In other words, we 
choose uji > and k x i > 0. Moreover, for definiteness we 



fix k w < 0. As a consequence, (38) and (39 1 are satisfied 
for k u > —k w . As to w/, it satisfies the asymptotic dis- 



persion relation u^uif — \s. 2 c 2 . A similar analysis can be 



performed for the solution which is regular at u = u+, 
which we denote by reg . Similarly, we denote by F^ 

and by F^ reg the singular and, respectively, regular so- 
lutions at 

Introducing the angle of emission 9, defined by k x — 
|k/| cos(#), Eq. (41 1 can be written as 



0Jl_ 

' 2v 



-n cos(#) 



(42) 



In states 



Prior to the formation of the pulse the wave equation 



„2 *! 



(43) 



whose independent positive frequency plane wave solu- 
tions are trivially given by (with uji > 0), 



exp (zk ; • X; - iuiti) 



(44) 



where the obvious dispersion relation holds: 



and where k\ = k 2 



(45) 



C. Quantization 

We can separate the monochromatic solutions 
{F + , F+ eg \ corresponding to all possible values of k;, into 
positive and negative frequency components {/ ki ,/k,}- 
These form a complete set of solutions which are mutu- 
ally orthogonal w.r.t. the index k; relative to the Klein- 
Gordon product in the retarded variables: 



(*,$) 



dudydz I 4" 



,5$ 
dw 



dw 



(46) 



The above product is a function of w (due to the factor 
n 2 (xi — vti) it is indeed impossible to define in a standard 
way a conserved inner product for Eq. Q). Expand- 
ing an arbitrary solution $ of Eq. ([!]) over the modes 
{/ k , /£ } we quantize the field by promoting $ to a field 
operator according to the expansion 



<$> = J d 3 kt (a ki / ki +</£,) 



(47) 



where the creation and annihilation operators a k; , a£ 
satisfy the usual commutation relations 



(48) 



Comparison between the creation and annihilation oper- 
ators of the in and out states will allow us to calculate 
the Bogoliubov coefficients relating one set to the other, 
thus leading us to evaluate the average number of emitted 
quanta. 



D. Greybody factor 

In our four dimensional problem, a priori we cannot 
neglect the effect of backscattering, i.e., the fact that, 
once emitted near the horizon, photons can undergo with 
a certain probability, reflection back in to the horizon 
due to the presence of a non-vanishing potential that 
they encounter in their propagation [2]. This originates 
the so-called greybody factor, which can be calculated 
as the square modulus of the transmission coefficient. 

Consider a massless scalar field propagating in a dielec- 
tric medium and satisfying the massless Klein-Gordon 
equation in the pulse metric (7 1, and let us call g(x) the 
component g TT appearing in ( 8 ) . We get 



n 2 (x) 1 
c 2 g(x 



d 2 $-n(x)d x 



-^g(x)d x <S>)-d 2 <S>-d 2 <S> = 0. 
n[x) J y 

(49) 



7 



We consider solutions of the form 

$(t, x, y, z) = e iuT tp(x)e- ik « y - ik * z . (50) 
As a consequence, we obtain an equation of the form 

{p{x)ip')' + cu 2 k{x)(p - q(x)cp = 0, (51) 
where the prime stays for the derivative w.r.t. x and 

n 



p(x) := -; k(x) := ~; q(x) := k\- 
n - ■ 
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eg 



The following change of variable: 

lk(u) 



s(x) 



p(u) 



du — , 
eg 



(52) 



(53) 



where s results to coincide with the tortoise coordinate, 
leads to the following Schrodinger-like form of the equa- 
tion for ip: 



d?<p 



where 



N , + (uj 2 - Q{a))tp = 0, 



Q{s) :=kl^. 

n 



(54) 



(55) 



It can be noticed that, as x — > x + , i.e. as s — > — oo, one 
obtains Q(s) —> 0, and that 



lim Q{s) = k 2 j_- 



1 — n 



(56) 



Q(s) turns out to be monotonically increasing from 
to the aforementioned constant value as s — > oo, and 
the transition from the zero value and the asymptotically 
constant value for a bump is very fast, in such a way that, 
as a plausible approximation for the calculation of the 
transmission coefficient, we replace the above potential 
Q(s) with a step-like effective potential with height Qca- 



d^ 



(57) 



In this approximation, we can refer to well-known results 
of standard quantum mechanics (see e.g. |34j . where the 
step barrier appears as a subcase of problem 37), in order 
to infer that: 

a) the transmission coefficient is zero if 



uj < 



n 



1 — n 



o„2 



b) the transmission coefficient is 



(58) 



(59) 



for 

w > V Qoo- 

We point out that the dispersion relation 



k sy , ~\~ k i 



(60) 



(61) 



implies that a) cannot occur. As a consequence, the 
Heaviside function which, in principle, should multiply 
the aforementioned |T| 2 , is always equal to 1. 

We have to translate the latter result in the laboratory 
frame; the simple substitution 



uj = {uii - vk x i)j, 



(62) 



where the label I indicates laboratory frame quantities, 
leads to the following expression for the greybody factor: 



r(w,,kj.) 

where 



4(uj t - vk x j)^(uji - vk x i) 2 - k\ql 
[(uji - vk xl ) + k±q ] 



(63) 



5o 



c 

n 



1 2 V 



One can also take into account that 

oji - vk xl =ui(l- n a ^- cos(9)j , 

so that 

r _ 4(l-n gcos(fl))|n g-cos(fl)| 

[l-n o ^cos(0) + K^-c< 

x 0h{ — (1 — no- cos(0)) — k±). 
qo c 



(64) 



(65) 



(66) 



Being 77 <^c 1, and v/c = (no+kn) 1 , fc e (0, 1), we obtain 

1 - n o -cos(0) - 1 - cos(0) + — cos(0), (67) 
c n Q 



n cos(0) ~ 1 — cos(0) , 

c n 



which leads to the conclusion that 

r ~ 1. 



(68) 



(69) 



It is worth pointing out that, on the grounds of the pre- 
vious result, a 90-degrees emission would not suffer any 
substantial suppression by backscattering. This is par- 
ticularly relevant in connection to recent measurements 
of 90-degree photon emission from RIP-induced horizon 



E. Thermal spectrum 



where 



In analogy with the corresponding black hole calcula- 
tions, in order to find the evaporation temperature of the 
RIP we must compute the expectation value of the num- 
ber operator of the outgoing photons evaluated in the in 
vacuum. In this connection, it is important to remark 
that, contrary to the standard black hole scenario, in the 
case of a RIP propagating in a nonlinear Kerr medium we 
have, strictly speaking, no real collapse situation. Nev- 
ertheless, we can "simulate a collapse" by building up 
a refractive index perturbation starting from an initial 
situation in which no signal is present in the dielectric. 

We perform a coordinate shift which moves the black 
hole horizon in u — 0. A photon starting at u < is 
trapped an cannot reach the front observer. As a con- 
sequence, states with u < cannot be available to the 
front observer and he will be lead to consider in his de- 
scription only the u > part of the F + solution, thus 
multiplying the latter by a Heaviside function Q(u). The 
Bogoliubov coefficients akjk,' an d Aqk,' relating, respec- 
tively, the positive and negative frequency components 
between the initial in and final out states are analysed 
in [Dj and satisfy the fundamental relation 



V 



k,k. 



(70) 



which is also proved in the appendix. Moreover, the ex- 
pectation value of the number operator of the outgoing 
photons in the initial vacuum state is 



(0 in\N% t \0in)=J2Ww 



(71) 



T 



1 



2nkBC 1 — -no cos I 



dn , 
du 



(74) 



By comparing this temperature with the temperature T + 
derived in the analogue Hawking model in the pulse frame 
(see Eq. (13l), we find the relation: 



T -n-^ c n co S (e) T+ - (75) 
For the same values of the parameters leading to 



Eq. (|2l|), for 9 = and for a typical value v ~ |c, we 



would obtain 



T ~ 78K, 



(76) 



which is again much greater than the values of T for typ- 
ical acoustic black holes. Moreover, in the specific case of 
a shock front, this temperature may increase significantly 
so that T ~ 2000 K (cf. subsection [TTC I . 



Equation. ( 75 1 gives the correct transformation law for 



the temperature in going over from the pulse frame to the 
laboratory frame. Indeed, to find how the temperature 
transforms, start from Wien's displacement law which 
gives the wavelength of maximum emission of a black 
body as a function of the temperature 



A ma „T = b = 2.9 x 10~ 3 m x K. 



(77) 



Converting to the frequency of maximum emission gives 



b 



(78) 



Then, by using the completeness relation for the Bogoli- 
ubov coefficients and taking into account backscattering: 

£(Kvl 2 -l/WI 2 ) =r(^,k x ), (72) 

k ; ' 



and Eq. (70), we obtain the following thermal-like dis- 



tribution written in terms of asymptotic (physical) fre- 
quencies: 



exp 



( hui \ 

\k b T ) 



(73) 



Now, under the boost connecting the lab frame to the 
pulse frame the frequency transforms according to the 
relativistic Doppler formula in a medium with refractive 
index uq 



lo = ujij 1 no cos I 

V c 



(79) 



where lo is the frequency in the pulse frame, and 9 is the 
emission angle with respect to the direction of motion 
of the pulse in the lab frame. Combining Eq. (79) with 
Eq. (1781) gives Eq. (1751). 



IV. EFFECTS OF OPTICAL DISPERSION 

In the framework of analogue models for black holes 
(dumb holes, bee models) there is a number of studies 
devoted to the analysis of the actual generation of 
Hawking radiation (see e.g. |33 [35H3H1 SOHII]). In 



particular, a suitable dispersion law, which involves the 
fourth power of the momentum k in the comoving frame 
is considered, since the original model by Unruh [35] - 
Calculations, in particular in presence of a black hole - 
white hole system, are quite involved and require also 
numerical simulations. 
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Black Hole - Pulse Frame 


Quantum Field Theory - Lab Frame 


propagation equation (eikonal approximation) 


propagation equation (eikonal approximation) 


1 


1 


appearance of a bh horizon 


onset of the surface s.t. 1 — n(u) - = 


1 


1 


logarithmic singularity in the phase of quantum modes 


logarithmic singularity in the phase of quantum modes 


I 


I 


thermal spectrum with temperature T+ 


thermal spectrum with temperature T 



TABLE I: Comparison between the analogue Hawking and QFT descriptions of quantum vacuum emission by a moving RIP. 



We do not deal herein quantitatively with the problem of 
mode conversion, and we do not consider the problem of 
black hole lasers (see [51 l3"3l |3"5TI4"3"] ). which will be taken 
into account in further developments of the present 
model. Still, we aim to point out the main physical 
consequences due to optical dispersion, and we shall 
stress that optical dispersion behaves as a fundamental 
ingredient for dielectric models. In this respect, there is 
agreement with dispersive models in sonic black holes. 

We start by recalling how optical dispersion affects the 
dispersion relation for the quantum electromagnetic field 
in linear homogeneous dielectric media (see e.g. [44] ) : 



n\uji) uf - k ; 2 c 2 = 



(80) 



This dispersion relation is obeyed by the monochromatic 
components of the field, and in [H] phenomenological 
quantization of the electromagnetic field in a dielectric 
medium is justified on the grounds of a more rigorous 
approach. In the presence of optical dispersion 



/>oo 

$(xj,tj)=/ duji<j)(uji,xi,ti), 
Jo 



(81) 



where (f>(u) h x h ti) = <p+(wi, xi)e- iwitl + <p-(w h xi)e iuit ' 
are the monochromatic components. Then it is easy to 
show that, even by allowing a spatial dependence for the 
refractive index 



V<p(wi , Xj ) +u t 



p(w l) x,) = 0. (82) 



Moreover, in the case of an homogeneous medium, by 
letting (p(u),x) cx exp(ik/ • x), one finds that eqn. (82) is 
satisfied only if 



By keeping into account a so-called 1-resonance model, 
we use the Sellmeier equation 

, ,2 

,2/ 



n\uJi) = 1 + 



(84) 



where ujq is the resonance frequency and u c is the cou- 
pling constant, which is also called plasma frequency. It is 
worth noting that the function n 2 (w;) is a Lorentz scalar 
and it is invariant under boosts. (84) amounts to a quar- 
tic equation in w/: 



where w 2 



K+ w 2 k ; V=0, (85) 
j 2 is the longitudinal resonance fre- 



quency. Solutions are 
1 



2„2 



k ; 2 c 2) 



4 W 2 



k;V 



(86) 

Then one finds two branches, which occur for < u>i- < 
lu and ui i < u>i + < oo respectively, with a forbidden 
region in the interval (ojq,u>l)- 

It is worth pointing out that, even if more complete ex- 
pressions can be provided for the dispersion relation, for 
practical purposes it is very useful the so called Cauchy 
formula 



n(u>i) =n + B u>?, 



(87) 

where hq is the would-be refractive index in absence of 
optical dispersion and Bq is a suitable constant; this ap- 
proximation works well for fused silica in the visible fre- 
quency interval, and can also be obtained from the more 



complete formula (84) in the formal limit as ui — > (i.e., 
from a physical point of view, as lo\ <C uiq). For our aims, 
a posteriori, its limited bandwidth validity is not a real 
problem, because it turns out that also Hawking effect 
takes place in a limited frequency window. 



■ u>i- 



n(u>i) 



(83) 



A. Dispersion and RIP 



i.e. only if each monochromatic component travels at its 
phase velocity (as obvious). 



We note that the dispersion relation for a linear homo- 
geneous dielectric medium is quite different from the one 
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which has been hitherto considered in the sonic model by 
Unruh and its subsequent developments (see also analo- 
gous comments in Ref. 05]). 

We may also wonder if the above relation can be trusted 
in our case, as we started from a nonlinear medium which 
is also non-homogeneous because of the RIP. We recall 
that we assume to be involved with linearized quantum 
fields around our effective geometry, and then it is rea- 
sonable to explore what happens in presence of optical 
dispersion by considering linear dispersion effects. More- 
over, we make the following ansatz: 



n(xi - vt u ui) = n (ui) + n f(x t - vti), 



i.e. we neglect optical dispersion effects in the RIP, and 
keep trace of it only in the background value n 0l which 
does not depend on spacetime variables, and passes from 
the status of constant to that of a function depending on 

In order to justify a posteriori this ansatz we can take 
into account what happens in the usual modelization of 
dielectric media when inhomogeneities are considered. A 
possibility is to account for inhomogeneities by includ- 
ing space-dependent density, polarizability, resonant fre- 
quencies, and so on (see e.g. [H]). By neglecting dissipa- 
tion in a polariton model, one can obtain a dielectric sus- 
ceptibility which depends on space and frequency, in such 
a way that the Sellmeier formula changes only because 
of an explicit dependence of lo c ,cjq on space variables. 
In our case, because of the properties of fused silica in 
the visible region, where the Cauchy formula works well, 
we can even neglect the spatial dependence of loq (which 
should be considered much greater than wj and uj c ), and 
then, recalling that u :— xi ~ vti, on the grounds of Ref. 
[4"6] we assume that 

u, 2 c (u)=uj 2 c0 +5oj 2 c1 (u), (89) 

where S <C 1 is a small parameter; if we also define 



n {uji) 



(90) 



the Sellmeier equation ( 84 1 then leads to 



n(ui,u) 




(91) 



where we have used the Cauchy approximation and ne- 
glected terms order of Surf. This validates our ansatz 



(88), with straightforward identifications. 

Another possibility is to adopt in the visible frequency 

region the following perturbative ansatz: 



n + nA x (u) + (B + nB^u))^ 

n + B u? +r ) A 1 (u), (92) 




2 4 

k(10 8 l/m) 

FIG. 2: (in color online) Branches relative to the dispersion 
relation ( 96 1 in the comoving frame. The dashed line corre- 
sponds to lo = — vk x . 



where no + Bqlu 2 = no^) and where we neglect terms 
0(r]ujf). The latter approximation is very useful, because 
we can write 



n(uj h u) = n (uji) + rjf(u) = n(u) + B q uj? 



(93) 



which isolates the contribution of the optical dispersion. 
As for the dispersion relation, we obtain in the lab 



n(ui,u)(Jt = ±\ki\c, 



n(u)tui + B ujf ± \ki\c = 0. 



(94) 



(95) 



This can be carried into the comoving (pulse) frame, 
where one obtains 



n(x)j(w + vk x ) ± \/7 2 (fc x + -i 2 ^) 2 + k]_c 



-B al 6 {u + vk x ) 6 = 0, 



(96) 



i.e. the same dispersion relation as in absence of optical 
dispersion except for the cubic term oc B . We shall 
consider in what follows the 2D reduction, where k± = 0. 
Graphical solutions of the above equation are displayed 
in figure [2] Note that 



oji = y(w + vk x ), 



(97) 



and then positive frequencies in the lab correspond to the 
region above the dashed line 



CO 



-vk x . 



(98) 



B. Dispersion and horizons 



The definition of horizon which is suitable for this kind 
of situation, where one cannot recover a metric by rein- 
terpreting the dispersion relation as an eikonal approxi- 
mation of the wave equation (due to the presence of more 
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than quadratic terms) , can a priori involve phase velocity 
and/or group velocity; cf. e.g. a short discussion in [48] . 
Let us first consider what happens when one involves 
the phase velocity, indicated as in the lab. reference 
frame and as v v in the comoving frame: 



V v u 



0, 



(99) 



In the lab. frame the calculation is elementary, and leads 
to 



v, 



±- 



~n(uji,u) ' 



which, due to ( 99 ) , implies the horizon condition 

c 



n(u>i,u) 



v 



(100) 



(101) 



it is worth noticing that condition ( 101 ) is a straightfor 



ward generalization to the dispersive case of the condition 
obtained in the absence of dispersion. It is also impor- 
tant to point out that the horizon condition in the lab in 
presence of dispersion can be also expressed as follows: 



1.457 



1.453 



1.449 




0.8 0.9 1 

wavelength (fxm) 



1.1 



FIG. 3: (in color online) Prediction of the Hawking emis- 
sion spectral range accounting for the RIP velocity and the 
material refractive index n. The two curves show, for the 
case of fused silica, no (in blue) and no + Sn (in green) with 
5n = 1 X 10" 3 . The shaded area indicates the spectral emis- 
sion region predicted for a Bessel filament that has an effective 
refractive index n = c/vb where the Bessel peak velocity is 
vb =Vq/ cos 8 with cone angle 6 = 7 deg and vq = du/dk is 
the usual group velocity. 



1 - w)- = 0, 
c 

which e.g. in the Gaussian model provides 



u± = ± 



n (oJi) 



V 



(102) 



(103) 



and, in general, a modified horizon condition appears: 



n (u!i) < - < n (uji) + r). 



(104) 



This condition is displayed in Fig. [3] (where wavelength 
replaces uji) and the velocity of the RIP is chosen to 
match the experimental conditions of Ref. [2B]. Explic- 



itly, by using the Cauchy formula (87), we obtain (we 



consider only positive frequencies) for rio+i] < ^, i.e. for 
a relatively small perturbation (see e.g. |26j ) 



1 S W i < y ■ 



i.e., in terms of wavelength 



.<A< 



n -r) 



(105) 



(106) 



If, instead, n + r\ > ^, i.e. for a substantially larger 
perturbation, we obtain 



< LJl < 

i.e., in terms of wavelength, 



VI- n o 



< A < oo 



(107) 



(108) 



For the parameters of Fig. [3] (see caption), Eq. (107) 
is verified for rj > 0.013, a rather large yet not com- 
pletely unrealistic refractive index variation. We may 
evaluate the influence of material dispersion in relation 
to the event horizon condition. In the absence of dis- 
persion a horizon is created only if v is tuned with ex- 
treme care such that Eq. ( 104 ) with dispersionless n is 



satisfied. Bearing in mind the typically small values of 
Sn - 10" 3 - 10" 4 (e.g. in fused silica, n 2 ~ 3 • 10" 16 
cm 2 /W and / ~ 10 13 W/cm 2 ), this would be no minor 
feat. 

Note that the Hawking photons will therefore in gen- 
eral be emitted only in a bounded spectral window. This 
is somewhat different from the dispersion-less case in 
which, once v is properly tuned so as to achieve the hori- 
zon condition, all frequencies are simultaneously excited, 
and, a posteriori, it can also justify why one can limit 
oneself to adopt the Cauchy formula instead of a more 
complete Scllmcicr. The aforementioned finite spectral 
window for photon emission in turn implies that in the 
dispersion-less case one should expect to observe the com- 
plete black-body spectrum predicted by Hawking. Con- 
versely, in the presence of dispersion in our analogue 
model, only a limited spectral region is excited and the 
black-body spectral shape will not be discernible. More- 
over, one has to take into account that, even in homo- 
geneous transparent dielectrics, optical dispersion affects 
the spectral density of photons which is actually to be 
taken into account, by introducing a phase space fac- 
tor multiplying the standard Planckian distribution term 
which depends on both the refractive index n(u>) and on 
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the group velocity v g (u>) [77] : 

/ Huj 3 \ n 2 (u) 



p(u>) = 



exp 



1 



(109) 



This formula is intended to hold true in the reference of 
the thermal bath. These optical dispersion contributions 
amounts to a sort of greybody factor arising from the 
interaction of photons with the dielectric material. See 
e.g. [37] • This is per se a sufficient reason for expecting 
deviations from the standard Planckian distribution. A 
very naive inclusion of dispersion in our case would also 
lead to T = T(oj), which would make even less plausible 
a pure Planckian spectrum. 

Regarding the width of the spectral emission window 
we note that this is determined by the value of n. For a 
typical case we consider a Gaussian pump pulse in fused 
silica, n Q ( 800 n m) = 1.467 and we find that for Sn = 



10~ 3 , Eq. (|104|) is satisfied over a ~15 nm bandwidth at 
435 nm. This window may become substantially large in 
lower dispersion regions as shown in Fig. [3] 

The phase velocity horizon condition may also be cal- 
culated in the comoving frame, where one obtains v v = 
if u = (for k x ^ 0), i.e. for v > 



n (y((J + Vk x ),x) ^(fc^a^O.fc^O} = 



(110) 



Notice that w = co(k x ,x) is solution of the 2D reduction 
of eqn. ( [96] ). The above condition requires that, as in 
absence of dispersion, no < f ■ Moreover, if also no + f? < 
-, then two disconnected regions are obtained: 



Vv 



n 



- n - rj 



< k x < 



V 



where only the latter region corresponds to U[ > 0. These 
are shown in Fig. [4] Of course, the latter condition is 
equivalent to condition ( 105 1 in the lab frame. 

If, instead, no 4 
is obtained, with 



' k x < 



n 



'Bo 



(111) 



> 



then a unique connected region 



~ n o 



< k x < 



(112) 



For a further 



This condition is consistent with ( 107 1 
comment, see below. 

One could also a priori consider a different definition of 
horizon, involving the group velocity instead of the phase 
velocity: 



vi.. 



V v a 



0. 



(113) 



where the first equality holds in the lab and the second 
in the pulse frame. For simplicity we consider a 2D re- 
duction of our model, so that k± = 0. One obtains 



v ^ 



dn(uji,u ) 
duj t 



-OJl. 



(114) 



-4 -2 2 4 

k(10 7 l/m) 

FIG. 4: (in color online) The dispersion relation (96) in the 
comoving frame for a 2D model at a fixed point x (only the 
branch relevant to horizon formation is shown). The shaded 
regions delimit the bands representing the values of k x for 
which a phase horizon exists. See Eqs. (110-111). 



In the case of the Cauchy formula, we find 

- = n(ui,u) + 2B ujf = n g (uji,u), 



(115) 



where n g indicates the group refractive index. In this 
case, in place of ( 104 1 , one would obtain 



n g {uji) < - < n g (uji) + rj, 



(116) 



which presents a clear shift with respect to the window 



corresponding to (104 1. Explicitly, by referring again to 
([87], we find 



V 7 ! -no-V \/f -no 

S < — == . 



i.e., in terms of wavelength 
2ircV3Bo~ 



< A < 



n 



n - n 



(117) 



(118) 



It is easy to realize that coexistence of phase and group 
horizons is allowed only if 



fBo 



i.e. if 



> 



(119) 



(120) 



Note that experiments, at least if the latter condition 
is not satisfied, should be able to distinguish between 
the two horizon conditions for phase velocity and group 
velocity given above. For example, in the experimental 
conditions of Ref. [26] used also in Fig. |3j condition ( 120 ) 
implies n > 0.009. However the maximum n obtained was 
only 0.001 and indeed the experiments clearly reveal 
emission in correspondence to the phase horizon alone. 
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It is possible to show that, if rj > ^ — no, there is 
a region in the (w, k Xl a;) -space where no group horizon 
appears because of the lack of any real zero in the deriva- 
tive of the dispersion equation. In order to realize this 
fact analytically, we can proceed as follows. Let us con- 
sider in the comoving frame the dispersion equation in 
the following form: 



D = D+D- = 0, 



(121) 



where 



D 



± - 



n(x)j(u!+vk x )+B ^ 3 (uj+vk x ) 3 zLj(k x +^ uAc = 0. 

c z 

(122) 

The dispersion relation uj{k x ,x) solves D = 0. The only 
dispersion equation which can lead to a group velocity 
horizon is the solution w_ (k x , x) of D_ = 0. Instead, the 
other branch u>+(k x , x) which solves D + = is monoton- 
ically decreasing. We find 



dcj(k x ,x) _ dD (dD\ 1 
dk x dk x \ dw 



(123) 



As a consequence, a group velocity horizon can emerge 



only by solving the system D— = and 



easy to show that = leads to 



dD- 



0. It is 



dk\ r 



k:\. 



n(x). (124) 



It is immediate to see that, if tiq + r\ > £ , as in the 
case where a unique connected r egion is found for t he 



phase velocity horizons (cf. eqn. (112)), k x ^ in (124) is 



complex valued for all x such that n(x) > ~, for which 
then group velocity horizons disappear. 

We may also qualitatively describe the behavior of 
a monochromatic or quasi-monochromatic wave near a 
phase horizon. In the comoving frame, an horizon is ap- 
proached only by waves that travel with v v > 0. Let us 
consider a wave traveling towards the white hole horizon; 
as it reaches the RIP, the refractive index increases, and 
then it slows down. This effect is enhanced by dispersion, 
because it implies a further increase of the refractive in- 
dex. As the wave is as near as possible to the phase (white 
hole) horizon, it nearly stops, but, as remarked by e.g. 
[36j . it cannot stop indefinitely, but, rather suffers reflec- 
tion (mode conversion). This happens for all frequencies 
belonging to the allowed windows indicated above. The 
previous analysis is corroborated by the study of null 
geodesies, in the presence of optical dispersion, which 
is carried out in |49j . Indeed, therein it is shown (in a 
more general geometric setting) that only geodesies with 
Vu, > undergo a slowing down process but, instead of 
suffering a process of trapping, as happens in absence of 
optical dispersion, they are reflected away from the trap- 
ping point. 

In the case of a wave packet, the dynamics appear to be 
more complex, but, due to its being a superposition of 



monochromatic components, we can infer that what hap- 
pens should be a sort of 'remastering' of the wave packet 
by the horizon in the following sense: frequencies beyond 
the allowed window are not affected by the presence of the 
horizon, whereas frequencies in the aforementioned win- 
dow are selectively 'bounced back' by a mode conversion 
mechanism. In this sense, we can appreciate the peculiar 
behavior of the phase horizon as a sort of semipermeable 
membrane, or, even better, as a sort of mirror with a 
selective reflection bandwidth, which transmits only the 
frequencies not belonging to the given windows. The re- 
flected frequencies will superimpose maybe even in the 
form of a wave packet. What is expected is that, in gen- 
eral, an incoming wave packet is converted into a reflected 
and possibly very broad and irregular wave train. 
The action of a group velocity horizon appears to be dif- 
ferent, in the sense that, although the existence of a win- 
dow of allowed frequencies seems to be analogous to the 
one of a phase velocity horizon, we have to point out 
that in the present case u>i refers to the carrier wave in 
the wave train. The group velocity horizon appears to be 
less selective than the phase velocity one: it limits itself 
to cause a reflection of wave packets which have carrier 
frequencies in the allowed window, without distinguish- 
ing between monochromatic components composing the 
packet itself. The outgoing wave train is expected to be 
still in the form of a compact and relatively undisturbed 
wave packet. Both numerical simulations and experimen- 
tal results could be able to reveal this different character 
of the two above horizon versions, and in particular mea- 
surements should be able to discern which definition is 
relevant for the physics at hand. It is worth pointing out 
that, in Ref. [26], the main role in photon production ap- 
pears to be related to the phase velocity horizon rather 
than to the group velocity horizon. 



V. CONCLUSIONS 

We have studied the photon production induced in 
a dielectric medium by a refractive index perturbation, 
which is created in a nonlinear medium by a laser pulse 
through the Kerr effect. The pulse has been assumed to 
have constant velocity v. In the pulse frame, we have 
investigated the analogous metric, and we have shown 
that two Killing horizons appear (black hole and white 
hole) for a generic but static dielectric perturbation. 
Despite the fact that the analogue metric is determined 
up to an overall conformal factor, the temperature is 
conformally invariant. Then we have taken into account 
quantum field theory in the lab frame, in a way that 
is independent of the analogous geometric framework 
we have again shown that there is a photon production 
with a thermal spectrum, corresponding to the Hawking 
effect. These findings are summarized in Table [I] 
We have also provided the transformation law of the 
temperature between the given frames and taken into 
account the effects of dispersion on the horizon condition. 
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Appendix A: Formal mapping to an acoustic black 
hole 



It is interesting to notice that our black hole metric in 
the pulse frame can be mapped into a form of acoustic 
black hole metric, provided that suitable identifications 
arc made. In particular, wc refer to the acoustic black 
hole metrics taken into consideration in |50| . and limit 
ourselves to the analysis of the x — t-part of the metric 
(indeed, only 2D metrics are studied in [50 ). In par- 
ticular, we are looking for a transformation allowing to 
carry 



ds 2 

ab (2) 



{<? - n 2 v 2 )dt 



7 



2v(l - n 2 )dtdx + n 



dx' 



(Al) 



\ Ml 



r/ 



FIG. 5: Penrose diagram for the analogue metric jlj. 



of •")() Fig. 28]. It is also straightforward to verify that 
the surface gravity can be calculated also by means of 
the formula 



d 

di {c - v) 



(A9) 



into the form 



dsi 



fl 2 [-{c 2 - v 2 )dt 2 - IMtdx + dx 2 ] 



where v plays the role of fluid velocity (in general depend- 
ing on t, x) and c is the local speed of sound, assumed 
to be constant as in [5U]. This re-mapping of the op- 
tical metric into an acoustic one is implemented by the 
identifications 



v = 

n 2 = 

dx 
dx 



7 v- 
1 



n 
1 



l 2 n 2 -^ 



"7 



(A3) 
(A4) 

(A5) 
(A6) 



It is straightforward to check that the horizon condition 
v = c is equivalent to the condition n = ^, i.e. the 
solution of ([9]). By defining 



du = dt — 



dw = dt 



dx 
c + v 
dx 



(A7) 
(A8) 



as in [SO], one can also easily draw the Penrose diagram 
of our spacetime (see Fig.[5|, which is analogous to that 



(cf. m\)- 



Appendix B: the Hawking temperature and its 
conformal invariance 



We consider the Euclidean version of the metric and 
find that near the horizon the ct — x part of the metric 
behaves like a cone which becomes a flat plane only if a 
special choice of an angle parameter, to be related to the 
inverse of the Hawking temperature, is chosen. 

The following rescaling to dimensionless variables is 
performed: 

2?r 

ct G (0,/3) h-> V = -j ct G (0,2tt) (Bl) 

2tt / 2tt \ _ . 

x G (x + ,oo) i-> x = —x G I — a;+,oo I , (B2) 

where r stays for the Euclidean time. Moreover, the 
following (local) diffeomorphism on the ct — x part of the 
metric is introduced: 



d,S(2) — 



2 r 



1 



n 2 {x) 



g TT {x)d(c 2 



1 



9tt{X) 



-dx' 



H 2 (y) (y 2 d^ 2 + dy 2 ), 



(B3) 
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whence 



1 



9tt(x) 
2 



n 2 {x) 
1 (dx 
g TT (x) \dy 



= H 2 (y). 



(B4) 
(B5) 



The above diffeomorphism is such that the ct — x part of 
the metric near the horizon can, at least locally, be made 
conformal to a 2D plane with a suitable choice of j3, with 

conformal factor (S^j H 2 (y). As a consequence, one 
finds (taking the positive sign) 



2tt dx 
Jdy 



1 1 



y n(x)' 



tt{x), 



i.e. 



y = Aexp 



2tt 

J 



dzn(z) 



9rr(z) 



(B6) 



(B7) 



where A is an integration constant. The diffeomorphism 
is defined to be regular if we can include in the manifold 
also the point y = 0, and then we have to require 

lim H 2 (y)= hm \g TT (x(y)) = h 2 , (B8) 

where h\ is a finite positive constant. The above limit is 
equivalent to 



lim 



1 



x-*x+ y 2 {x) 
Near the horizon one finds 



9rr{x). 



y 2 ~(x-x + y^^^^ 



and, as a consequence, one has to choose 



c 2 1 



1 



v 2 7 2 



=: he fih 



(B9) 



(BIO) 



(Bll) 



The above method can be used also to confirm that 
the temperature does not depend on the (static) confor- 
mal factor. Indeed, an overall conformal factor tt 2 (x) 
which is finite and non- vanishing at the horiz ons ( 
lim x _j. x± fl 2 (x) — b 2 , > ) does not modify Eq. (BIO I, 
as may be realized using Eqs. (B8| and (B9|. 



Appendix C: Exact analysis of equation ( 26 ) near 
the singularity 

We confine ourselves to discuss the root u+, since all 
considerations can be trivially extended to u_ . Then, we 
rewrite Eq. ( 26 ) as 



A"(u) + J^La'(u) + ^^Afu) = 0, (CI) 
u — u+ u — u+ 



where we have introduced the functions P(u) and Q(u) 
(holomorphic in the disk \u — u + \ < \u + — 



(C2) 
(C3) 




In a neighbourhood of w + Eq. ( 26 ) writes, to leading 
order in u — u + , as 



A"(u) 



1 



vn'{v,-\ 



-A'(u)+ 



(k 



2vn'(u-\-)(u — 



-A{u) = 0, (C4) 



where the prime denotes derivation with respect to the 
lab variable u. The indicial equation is a(a — 1) — 

1 



2ik„ 



a — = 0, with roots ot\ = and 012 



Therefore, in the neighbourhood of m + , Eq. (26 1 has two 
linearly independent solutions of the form 



vn f (u+) ' 



A^\u) = (u- 



n=0 



(C5) 



i = 1,2, where the series define holomorphic functions 
in the disc \u — u+| < \u + — U—\ and whose coefficients 
can be obtained recursively from the equation c~n (oii + 

n)(a i +n-l+p )+Y^Zl{{o L i+ r )Pn-r+qn~r-l)Cr ) = 0, 

where the Pk, Uf. (k — 0,1,2,.. .) are the coefficients of 
the expansion about u + of P(u) and Q(u) respectively. 
Of particular interest is the solution ( 25 ) corresponding 
to a — a-2, which for u > u + has the form 



F a2 (u,w,y,z) = ^(^e^'^+i 



log(u—u + )-\-ik w w-\-ik y y-\-ik z z 

(C6) 

where £ (u) is holomorphic in the neighbourhood of u = 
u+ and vanishes as u — > u+: indeed, we have reabsorbed 
in it the factor (u — u + ) associated with the real term in 
a2, in such a way that 



£(u) = (u-u + )r)(u), 



(C7) 



where ry(u) is holomorphic in the neighbourhood of u = 
u+ and i](u) = Co + c\(u — u+) + . . . as u — > u+. We 
stress the presence of the logarithmic divergence of the 
phase as u approaches u+ even in this approach. Also, 
by comparison with the study of the solutions contained 
in section III A, we can infer that the above solution cor- 
responds to the (exact) expansion of F^(u,w,y,z) near 
u = u+. 



Appendix D: Analytic continuations 

We are interested in the Bogoliubov coefficient o^', 
which relates the positive frequencies between the initial 
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in state: 



pin ik x x-\-ik y y-\-ik z z — icjt 



(Dl) 



and final out state F^: 



F+ = 9{u - u+ )^{u)e l[ayo&[u - u+)+k ™ w+k y y+k * z \ (D2) 

where we have introduced the shorthand notation a = 
,, 2c x k,„ and, in agreement with the results obtained in 
Appendix [Cj we have introduced also the analytic part 
It will be evident that, at least in the large fre- 
quency limit, this cannot affect the thermal charac- 
ter of particle emission. 

Let us start with the computation of , which, apart 
from a factor which will not affect our goal, which consists 
in the deduction of (D7|, is given by 



/>oo 

a kk , cx<5 2 (k ± -kl) / ^{u)u ia t 
Jo 



du, 



(D3) 



where we have shifted the variable u so that u + is mapped 
on 0. Since 2k' u = (k' x +uj' /v) > 0, we see that the second 
exponential factor is rapidly decreasing at infinity when 
lm(u) < 0. Consider the path T starting from to R > 
along the real line, then following the arc of radius R 
clockwise until — iR, and finally coming back from — iR to 
along the imaginary axis. As the integrand is analytic 
inside the region bounded by the path, the integral along 
r vanishes for any positive value of R. Then, taking 
the limit R — > +oo we then see that the integral along 



the positive real axis is equal to the integral along the 
negative imaginary axis and 



«kk' oc e 2 



^ 2 (k ± -kl) 



-it 



(D4) 



We can consider the limit as k' u ^> 1 and approximate £ 
for small values of its argument. Then we obtain 



a k k' oc e 



^ 2 (k^-k- 



(K) 



2 + ia 



T(2 + ia), 



(D5) 



where T is the Euler Gamma function. 

On the other hand, in order to calculate /3 k k' j it is suf- 



ficient to revert the sign of k' u and k' in Eq. ( D3 ) . This 



time, the integral can be calculated along the positive 
imaginary axis, since now the term relative to the inte- 
gration on the arc vanishes for Im(u) > 0. Thus, rotating 
the path counter-clockwise, we obtain 



Ack' oc e 



(K) 



-T(2 + ia), 



(D6) 



By comparing Eqs. (D5) and (D6|, it is easy to verify 
that 



Okk' 



1/3: 



kk' 



(D7) 



which is Eq. (70) 
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